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Helimagnetic thin films: surface reconstruction, surface spin-waves and magnetization 
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Quantum properties of a helimagnetic thin him of simple cubic lattice with Heisenberg spin model 
are studied using the Green’s function method. We hnd that the spin conhguration across the him 
is strongly non uniform. Using the exactly determined spin conhguration we calculate the spin-wave 
spectrum and the layer magnetizations as functions of temperature T. We show the existence of 
surface-localized modes which strongly affect the surface magnetization. We also show that quantum 
huctuations cause interesting spin contractions at T = 0 and give rise to a cross-over between layer 
magnetizations at low T. 
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I. INTRODUCTION 

Recently, there has been a growing interest in magnetic 
properties of helimagnets, due to possible applications 
in spin transport properties usingmaterials at nanoscale 
such as thin hlms and multilayers [l| . In particular, inten¬ 
sive researches have been carried out to understand the 
role of skyrmions There is in addition a large num¬ 

ber of experiments which has recently been performed on 
thin hlms of helical magnets [iMl. 

The simplest model of the helimagnetic ordering is 
due to a competition between nearest-neighbor (NN) 
and next-nearest-neighbor (NNN) interactions, as dis¬ 
covered by Yoshimori 0 and Villain 0 : a spin in a 
chain turns an angle 0 with respect to its previous neigh¬ 
bor. There are many families of helimagnets due to var¬ 
ious kinds of interaction among them one can mention 
non collinear magnetic structures due to Dzyaloshinskii- 
Moriya interactions or to geometry frustration 0iI3. 
Low-temperature properties in helimagnets such as spin- 
waves and heat capacity have been exten¬ 

sively investigated. In spite of their long history, the 
nature of the phase transition in non collinear magnets 
such as stacked triangular XY and Heisenberg antiferro- 
magnets has been elucidated only recently |23l-[^. For 
reviews, the reader is referred to Ref. [26|. 

In this paper, we study a quantum Heisenberg helimag¬ 
netic thin film with the simple cubic (sc) lattice. The case 
of the body-centered cubic (bee) lattice has been recently 
studied [^. Surface effects in thin films have been in¬ 
tensively studied during the last three decades [H, [29| . 
However, due to complicated surface spin configurations, 
surface effects in helimagnets have only been recently 
studied: surface spin structures [s^, Monte Carlo (MC) 
simulations [3l| and a few experiments [32|, [33|. Heli¬ 
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cal magnets present potential applications in spintron- 
ics with predictions of spin-dependent electron transport 
in these magnetic materials |1Q|-[1^. This motivates the 
present work. 

We shall use the Green’s function (GF) method which 
has been initiated by Diep-The-Hung et al. for collinear 
surface spin configurations [sj. For non collinear mag¬ 
nets, the GF method has also been developed for bulk 
helimagnets 0 and for frustrated films [sHI, [36| . In he¬ 
limagnets, the angles between neighboring spins become 
strongly non uniform as seen below, making calculations 
harder. This explains the small number of microscopic 
calculations so far for helimagnetic films. 

The paper is organized as follows. In section II, the 
model is presented and classical ground state (GS) of the 
helimagnetic film is determined. We summarize there the 
principal steps used in the general GF method for non- 
uniform spin configurations. The GF results are shown in 
section III where the spin-wave spectrum, the zero-point 
spin contraction and the layer magnetizations are shown. 
Concluding remarks are given in section IV. 


II. MODEL, CLASSICAL GROUND STATE AND 
QUANTUM FORMULATION 

We consider a thin film of sc lattice of Nz layers, with 
two symmetrical surfaces perpendicular to the c-axis, for 
simplicity. The exchange Hamiltonian is given by 

= - X • Si (1) 

ihj) 

where Jij is the interaction between two quantum 
Heisenberg spins and Sj occupying the lattice sites 
i and j. 
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A. Surface spin reconstruction 

To generate a bulk helimagnetic structure, the simplest 
way is to take a ferromagnetic interaction between NNs 
Ji (> 0 ): and an antiferromagnetic interaction between 
NNNs J 2 < 0. If IJ 2 1 is smaller than a critical value I J 9 1 , 
the classical GS spin configuration is ferromagnetic [l 8 |- 
113 . Let us consider the case of a helimagnetic structure 
only in the c-direction perpendicular to the film surface. 
In such a case, we assume a non-zero J 2 only on the c- 
axis. This assumption simplifies formulas but does not 
change the physics of the problem since including the 
uniform helical angles in two other directions parallel to 
the surface will not introduce additional surface effects. 
The bulk quantum heliniagnets have been studied by the 
Green function method |^ . 

For the present model, the helical structure in the bulk 
is planar: spins are parallel in planes perpendicular to 
the c-axis and the angle between two NNs in the adja¬ 
cent planes is a constant and is given by cos a = 
for a sc lattice. The helical structure exists therefore 
if IJ 2 I > O. 25 J 1 , namely |J 2 |(bulk)= 0.25Ji. To calcu¬ 
late the classical GS surface spin configuration, we write 
down the expression of the energy of spins along the c- 
axis, starting from the surface: 

E = - Ji cos((9i - 6 > 2 ) - Ji [cos((92 - 6 >i) 

+ cos(6>2 - O3 )] + ... 

—J2 COs( 6 >i — O3) — J2 COs(^2 — O4) 

- J2[cos(6>3 - Oi ) + cos(6>3 - O3 )] + ... (2) 

where Oi denotes the angle of a spin in the i-th layer made 
with the Cartesian x axis of the layer. The interaction 
energy between two NN spins in the two adjacent layers 
i and j depends only on the difference ai = 9i — 

The GS configuration corresponds to the minimum of E. 
We have to solve by iteration the set of equations: 

BE 

— = 0 , for i = (3) 

BOLi 

The result is shown in Fig. [ 1 ] for Nz = 8 . Some re¬ 
marks are in order: i) the result is obtained by iteration 
with errors less than 10 “^ degrees, ii) strong angle vari¬ 
ations are observed near the surface with oscillation for 
strong IJ 2 I, hi) the angles at the film center are close to 
the bulk value a, meaning that the surface reconstruction 
affects just a few atomic layers (this is more clearly seen 
for thicker films not shown here). This bulk helical sta¬ 
bility has been experimentally observed in holmium films 
[ 3 ^ . An alternative method giving the same result is the 
numerical stee pes t descent method which is described in 
details in Ref. l35l . 



FIG. 1: (Color online) Angles ai . 07 in degree across the 

film for J 2 =- 0 . 6 , -0.5, -0.4, -0.35, -0.3 (from top) with Nz — 8 . 

spin coordinates defined as follows mill: the quanti¬ 
zation axis of spin Si is on its Q axis which lies in the 
plane, the r]i axis of Si is along the c-axis, and the ^i axis 
forms with rji and Q axes a direct trihedron (see Fig. [2]). 



FIG. 2: Local coordinates in a xy-plane perpendicular to the 
c-axis. Q denotes Oj — Oi. 

Expressing the Hamiltonian in the local coordinates, 
we obtain 

= - E 4iU(cos%-l)(5+5+ + 5-5-) 

<ij> I 

+ i(cos% + l)(5+5-+5-5+) 

+ i sineij {St + S-) i sineijS! {St + S-) 

+ cos%S'fS'|| (4) 


B. Analytical formulation 

To calculate physical quantities at finite temperatures, 
we shall use the GF method. To that end, we use the local 


Now, according to the theorem of Mermin and Wagner 
[HI continuous isotropic spin models such as XY and 
Heisenberg spins do not have long-range ordering at finite 
temperatures in two dimensions. Since our films have 
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small thickness, it is useful to add an anisotropic inter¬ 
action to stabilize the long-range ordering at finite tem¬ 
peratures. Let us use the following in-plane anisotropy 
between and S^: 

^a = - ^ (5) 

<hj> 

where Iij{> 0) is supposed to be positive, small com¬ 
pared to Ji, and limited to NNs. The full Hamiltonian is 
thus H = He Ha- The GS in the presence of Iij{> 0) 
can be determined in the same manner. Hereafter we 
take lij = Ii for any NN pair, except otherwise stated. 
It is only very slightly modified with the order of one or 
two degrees when Ii O.lJi. The small anisotropy does 
not therefore alter the main features shown in Fig. [TJ 
The general method has been recently described in de¬ 
tails in Ref. 113 . To save space, let us just briefly recall 
here the principal steps of calculation and give the results 
for the sc helimagnetic film only where they should be. 
We define the following two double-time Green’s func¬ 
tions in the real space: 

= «s+{ty,s-{t’)» 

= -i0{t-t')<[St{t),S-{t')]> (6) 

= «s-ity,s-it')» 

= < [S~{t),S~{t')] > (7) 

We need these two functions because the equation of mo¬ 
tion of the first function generates functions of the second 
type, and vice-versa. Writing the equations of motion 
of these functions and using the Tyablikov decoupling 
scheme to reduce the higher-order functions, for example 

« Sf,sntysr{t') »c^< >« s+itys-in » 

etc., we obtain the general equations for non collinear 
magnets p^ . 

We next introduce the following in-plane Fourier trans¬ 
forms: 


= ^ [ f dk, T dcoe-A^-^') 

J J B Z J — cxD 

Xflni.nj (ujAxv) (8) 

FkjitA) = X / / rfk, T f °° dcue-^A-*') 

J J B Z J — ( 


' BZ J —CO 

pikajy • (Flfc Rj) 


^ fnk,nj ^ 


(9) 


where uj is the spin-wave frequency, denotes the wave- 

vector parallel to xy planes and is the position of the 
spin at the site i. rij and rik are respectively the z- 
component indices of the layers where the sites R^, Hj 
and R/c belong to. The integral over is performed in 
the first Brillouin zone {BZ) whose surface is A in the 
xy reciprocal plane. For convenience, we denote rii = 1 
for all sites on the surface layer, rii = 2 for all sites of the 
second layer and so on. 

We finally obtain the following matrix equation 


M (cj) h = u. 


( 10 ) 


where M (cj) is a square matrix of dimension 
(21V. X 2N2:)^ h and u are the column matrices which 
are defined as follows 


/ 9l,n' \ 

flA ' 


h = 


9n,n' 

fn,n' 


,U = 


V fN^,n' / 


2{St)5en' \ 

0 


0 / 


where, taking h=l hereafter. 


(11) 


J 


( UJ 


0 

Bt 

ct 

Dt 

Et 

0 

0 

0 

0 

0 

0 

\ 


0 

UJ — Ai 

-ct 

-Bt 

-Et 

-Dt 

0 

0 

0 

0 

0 

0 




Dn 

Ey 

Bn 

Cn 

UJ + An 

0 

Bt 

Ct 

Dt 

Et 





-En 

-D- 

-Cn 

-Bn 

0 

UJ Afi 

-Ct 

-Bt 

-Et 

-Dt 




0 

0 

0 

0 

0 

0 



Bn^ 

Cn^ 

UJ + A^^ 

0 



0 

0 

0 

0 

0 

0 


~Bn^ 

~Cn^ 

~Bn^ 

0 

UJ — An^ 

/ 







where n = 1,2,. 

..,Nz, dn = h/Ji, and 

= -8J// < S', 

^ > (1 + dn 

- 7 ) 


Bt 

= 2Jf^ (S^) (cos6'„,„±i + 1) 

-2 < 

qZ 

^n+1 

> cos dn,n-|-l (dn 

+ ) 

Ct 

= 2Jf^ (S^) (cos6»„,„±i - 1) 

-2 < 

qZ 

^n-1 

^ COS Ofi^n— 

-1 (dn 

+ Jf^) 

Et 

= J 2 (S^) (cos6>„,„±2 - 1) 

1 

to 

< S ^+2 > cos6>„, 



Dt 

= J 2 (S^) (cos6>„,„±2 + 1) 

1 

to 

<st 

_2 > COsdn, 

n—2 

(13) 




where 
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Note that to use the above formulas, we have to apply 
the following rules: (i) if n = 1 then there are no n — 1 
and n — 2 terms in the matrix coefficients, (ii) if n = 2 
then there are no n — 2 terms, (iii) if n = then there 
are no n + 1 and n + 2 terms, (iv) if n = — 1 then 

there are no n + 2 terms. Besides, we have distinguished 
the in-plane NN interaction from the inter-plane NN 
one J^. 

III. RESULTS AND DISCUSSION 

Using the spectral theorem which relates the correla¬ 
tion function {S^S^) to the Green’s function [^, we 
have 

{S-S+) = lim^J J dKy j ^{9n,n'{i^ + is) 

- (14) 

where e is an infinitesimal positive constant and /3 = 
ks being the Boltzmann constant. Using the 
Green’s function presented above, we can calculate self- 
consistently various physical quantities as functions of 
temperature T. The magnetization {S^) of the n-th layer 
is given by 

(^n) = 

+ 00 

— oo 

(15) 

After some steps, we obtain 

= l + ( 16 ) 

^ i=l 

where n = and D 2 n-i(<^i) is the determinant 

obtained by replacing the (2n — l)-th column of M by 
u at As < 5'^ > depends on the magnetizations of 
the neighboring layers via cJi(i = l,...,2N;z), we should 
solve by iteration the equations ([T6|) written for all layers, 
namely for n = 1, ...,Nz, to obtain the magnetizations of 
layers 1, 2, 3, ..., Nz at a given temperature T. Note that 
by symmetry, < S'f >=< >, < >=< >, 

< S'! >=< >, and so on. Thus, only Nz/2 self- 

consistent layer magnetizations are to be calculated. 

The value of the spin in the layer n at T = 0 is calcu¬ 
lated by 


11 /*/* 

+^)(r = 0) = - + -y j dk,dkyY,D2n-l{uJi) (17) 


where the sum is performed over Nz negative values of 
uji (for positive values the Bose-Einstein factor is equal 
to 0 at T = 0). 

The transition temperature Tc can be calculated in a 
self-consistent manner by iteration, letting all < > 

tend to zero, namely cji —> 0. Expanding e^^* — 1 ^ /3cUJi 
on the right-hand side of Eq. m where + = {kBTc) ^, 
we have by putting (S'^) = 0 on the left-hand side, 

= 4 / / % ++ ( 18 ) 

There are Nz such equations using Eq. m with n = 
1 ,...,A; 2 . Since the layer magnetizations tend to zero 
at the transition temperature from different values, it is 
obvious that we have to look for a convergence of the 
solutions of the equations Eq. m to a single value of 
Tc. 

A. Results 

Let us take = J = 1 everywhere except on 

the surface where j{^ = Jg. We use d = hj/J for any 
NN pair, for simplicity. 

Numerically, we use a Brillouin zone of 100^ wave- 
vector values, and we use the obtained values (S^) at 
a given T as input for a neighboring T. At low T and up 
to ^ |Tc, only a few iterations suffice to get a conver¬ 
gence precision <1%. Near Tc, the convergence is much 
harder. We show below our results. 

We have calculated the spin-wave spectrum uj versus 
kx = ky for various values of J 2 in the case of a eight- 
layer film with an anisotropy d = 0.1. There are 8 pos¬ 
itive and 8 negative modes corresponding two opposite 
spin precessions. We can mention here the existence of 
acoustic surface modes which lie in the low energy region 
for Jg = 0.6 as seen in Eigl3] (middle) and optical surface 
branches which lie outside the bulk-mode energy region 
for Jg = 1.6 seen in Eigl3] (bottom), whereas no such 
modes exist in the case when Jg = 1 [Eig|3] (top)]. 

It is known that in antiferromagnets, quantum fluc¬ 
tuations give rise to a contraction of the spin length at 
zero temperature . We will see here that a spin under 
a stronger antiferromagnetic interaction has a stronger 
zero-point spin contraction. The spins near the surface 
serve for such a test. In the case of the film considered 
above, spins in the first and in the second layers have only 
one antiferromagnetic NNN while interior spins have two 
NNN, so the contraction at a given J 2 /T 1 is expected to 
be stronger for interior spins. This is verified with the 
results shown in Eig. |U When IJ 2 I/J 1 increases, namely 
the antiferromagnetic interaction becomes stronger, we 
observe stronger contractions. Note that the contraction 
tends to zero when the spin configuration becomes ferro¬ 
magnetic, namely J 2 tends to -0.25. 

We show the layer magnetizations in Eig. [5] in the case 
where J 2 /J = —0.7 and Nz = 8. Some remarks are in 
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FIG. 3: (Color online) Spin-wave spectrum versus k = kx — 
ky in the case where Nz = S and d = 0.1 for Js = 1 (top), 
Js = 0.6 (middle) and Js = 1.6 (bottom). 


order: 

(i) the shown result is obtained with a convergence of 
1%. For temperatures closer to the transition tempera¬ 
ture Tc, we have to lower the precision to a few percents 
which reduces the clarity because of their close values 
(not shown). 

(ii) the surface magnetization, which has a large value 
at T = 0 as seen in Fig. 01 crosses the interior layer mag¬ 
netizations at T 0.6 to become smaller than interior 
magnetizations at higher temperatures. This cross-over 
phenomenon is due to the competition between quan¬ 
tum fluctuations, which dominate low-T behavior, and 
the low-lying surface spin-wave modes which strongly di¬ 



J2 


FIG. 4: (Color online) Spin lengths at T = 0 for several values 
of J 2 with d = 0.1, Nz — 8: black circles, void squares, black 
squares and void circles are data for spins in first, second, 
third and fourth layers, respectively. 



T 

FIG. 5: (Color online) Layer magnetization as function of T 
for J 2 = —0.7 with d = 0.1, Nz = S: red circles, blue void cir¬ 
cles, magenta squares and black void squares are magnetiza¬ 
tions of the first, second, third and fourth layers, respectively. 


minish the surface magnetization at higher T. Note that 
the second-layer magnetization makes also a crossover 
at T 0.6. Similar cross-overs have been observed in 
quantum antiferromagnetic films 0 and quantum su¬ 
perlattices jdlj . 

Note that though the layer magnetizations are different 
at low temperatures, they will tend to zero at a unique 
transition temperature as seen below. The reason is that 
as long as an interior layer magnetization is not zero, it 
will act on the surface spins as an external field, prevent¬ 
ing them to become zero. 

Let us show in Fig. [6] another example of layer mag¬ 
netizations (without zoom at low T) up to temperatures 
close to the transition, for J 2 = —0.5. The convergence is 
rather good but it is difficult to get to Tc. We explain how 
to determine Tc by another way which is easier. As said 
earlier, each equation m for a given n gives a pseudo 
transition temperature Tcs as long as T is not close to 
the temperature where all layer magnetizations are very 
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FIG. 6: (Color online) Layer magnetizations as function of T 
for J 2 = —0.5 with d — 0.1,ds = 0.2, Nz = 8: red circles, 
green void triangles, blue triangles and magenta circles are 
magnetizations of the first, second, third and fourth layers, 
respectively. 



FIG. 7: (Color online) Example of how to determine the tran¬ 
sition temperature for J 2 = —0.5 with d = 0.1, ds = 0.2, 
Nz = 8 (see text): red, green, blue and magenta circles are Tcs 
determined from formula m for n = 1,2, 3, 4, respectively. 
Their intersection gives the transition temperature Tc — 2.33. 

small. To determine this temperature, we plot Tcs ob¬ 
tained at several temperatures. The convergence of these 
temperatures to a single one occurs when T = Tc. This 
is shown in Fig. [71 

B. Discussion 

Let us compare the results found in this paper for a 
thin film of sc lattice and those for a thin film of bee 
lattice studied in Ref. 

(i) both represent a strong non uniform spin recon¬ 
struction as a function of J 2 . Note that the critical value 
J 2 is —O. 25 J 1 in the sc case while it is — Ji in the bee 
case. So, the angle variation at the surface is not the 
same for a given value of J 2 in the two cases. 

(ii) both show a cross-over of layer magnetizations at 
low temperatures, however the order of the layer magne¬ 
tizations before as well as after the cross-over is not the 


same in the two cases. 

(iii) the zero-point spin contraction is different in two 
cases: the sc case shows the first-layer spin contracts less 
than the second, the second less than the third, the third 
less than the fourth (see Fig. 0]), while in the bee case 
the fourth layer contracts less than the others (see Fig. 4 
of Ref. 113 ). This is in agreement with the spin contrac¬ 
tions discussed in point (ii) above and can be understood 
by looking at the antiferromagnetic contribution to the 
local field at a spin of each layer: the smaller this con¬ 
tribution the smaller the contraction. Besides, the bee 
spins contract more strongly than the sc ones. 

(iv) the spin-wave spectrum is different in the two 
cases: in the case where surface interactions are the same 
as the bulk interactions, the sc spectrum does not have 
surface-localized spin wave while the bee spectrum has 
an acoustic surface branch ver y si milar to the antiferro¬ 
magnetic cases shown in Ref. \3^ . This is because the 
surface spins lack four NN while the sc spins lack only 
one NN. When surface interactions are smaller (larger) 
than the bulk ones the sc shows acoustic (optical) surface 
modes (see Fig. El). The bee case shows similar effects 
but at different values of J|. 

The above qualitative and quantitative similarities and 
differences are very important when one deals either the¬ 
oretically or experimentally with the films of different 
lattice symmetries. 


IV. CONCLUSION 

Surface effects in a helimagnet of simple cubic lattice 
with quantum Heisenberg spins have been investigated in 
this paper starting from the classical ground-state spin 
configuration which is exactly determined. The strong 
surface spin rearrangement is observed but it is insensi¬ 
tive to the film thickness in agreement with experiments 
performed on MnSi films [32| and holmium . We have 
calculated self-consistently physical quantities such as the 
spin-wave excitation, the spin length at T = 0 and the 
layer magnetizations as functions of temperature. We 
have shown that when varying the surface exchange in¬ 
teraction, we observe surface-localized acoustic and opti¬ 
cal modes which lie outside the propagating-magnon en¬ 
ergy band. These modes cause a strong deviation of the 
surface magnetization with respect to the interior ones. 
Another interesting phenomenon is the cross-over of layer 
magnetizations at low temperatures due to the competi¬ 
tion between quantum fluctuations and thermal effects. 
A comparison of the results found here with those for the 
bee case has been given. 
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